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We study multigap superconductivity, with strong angular variations of one of the gaps, as ap- 
propriate for certain iron-based high-temperature superconductors. We solve the gap equations of 
this model and find that the nodes or zeroes in the gap function present at T c - although purely ac- 
cidental - typically survive down to T = 0. Based on this result, we investigate the line of quantum 
transitions at which gap zeroes first appear. The peculiar "zero-point" critical scaling emanating 
from this line dominates quasiparticle thermodynamics and transport properties over much of the 
phase diagram, and supplants more familiar forms of scaling associated with the accidental nodes. 



I. INTRODUCTION 

The discovery of iron-based high-temperature super- 
conducting family 1 has altered the landscape of con- 
densed matter research. The underlying physics of iron- 
pnictides and the first high-T c family - the cuprates - ap- 
pear to be significantly different. Thus, after two decades 
of cuprate domination, the superconductivity research 
faces new major challenges and new important problems. 

The progress in our understanding of iron-pnictides 
has been swift 2 ^, but many important questions remain 
unanswered. In particular, the form of the order param- 
eter remains uncertain. There is growing consensus that 
the it belongs to the general class of the so-called s±, s' 
or extended s-wave stated, the basic dynamical origin of 
which can be understood from analytic renormalization 
group (RG) arguments^. This superconducting state is 
generated by the electron-electron interaction, originat- 
ing in Coulomb repulsion. This mandates a sign change 
between the gaps on the multiple sheets of the Fermi 
surface (FS) (interestingly, a similar order parameter ap- 
pears within the strong coupling approach to pnictidcs 7 ). 
In 122s - the best studied compounds of this family - the 
angle resolved photoemission spectroscopy (ARPES)£~— 
detects up to four relatively isotropic gaps, one for each of 
the two hole and two electron bands of the FS. This is in 
an apparent contradiction to the penetration dept h 11 ' 12 , 
thermal conductivity^ and specific heat dat a 14 ' 15 , which 
indicate gapless quasiparticle excitations - a natural in- 
terpretation is that these are associated with nodes in the 
gap somewhere along the (multiply-connected) FS. More- 
over, Refs. [TH and [TBI report such excitations both in the 
under- and over-doped phases of (Co,Ba)122, but they 



apparently disappear at optimal doping. This suggests 
a significant change in the gap structure as a function 
of doping. Similar doping evolution of the gap has been 
reported based on penetration depth measurements 1 ^. 

Such reading of the experiments receives some sup- 
port from theory. Random phase approximation (RPA)- 
based calculation o 17 ' 18 hint at a multitude of compet- 
ing nodal and nodeless states and possible transitions 
among them, tuned by material parameters like pnicto- 
gen height or impurity concentration. Furthermore, fluc- 
tuation exchange 1 ^ and numerical functional RG (fRG) 
studies 2 ^ - — also find significant angular variations and 
possible nodes in one or more of the gaps. A similar pos- 
sibility was discussed in a model of Ref. 23. In these 
studies, the nodal structure is typically induced by the 
strong orbital anisotropy of effective interactions within 
the extended s-wave state. Thus, such nodes are "acci- 
dental, " in the sense that they are not protected by any 
symmetry or topological considerations. Consequently, 
their presence or absence, and location on the FS are all 
affected by the interactions, temperature, impurity scat- 
tering, and the like. 

Both the numerical fRG and RPA studies, as well as 
Ref. |23| . have focused on the linearized gap equations 
and the ensuing pairing states at T — T c . However, 
since the gap equations form a complicated non-linear 
system there is no guarantee that the gap structure will 
remain unchanged, and that such accidental nodes or ze- 
roes, present at T = T c , survive as T — > 0. 

In this paper we consider a model of a two-band super- 
conductor, which has a uniform gap on one band (hole) 
and gap with angular variations on the other (electron) 
band. The second gap takes the form Aq + A' cos 4(9, 
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which obviously allows nodes for A' > Ao. As stressed 
above, both nodal and nodeless states belong to the na- 
tive A\ g representation of the lattice rotation group and 
the appearance of nodes is a direct consequence of the or- 
bital anisotropy of the effective interactions. This model 
emulates the physics of iron-pnictides and its different 
versions have been used in that context - — . A sum- 
mary of our results is as follows: i) The gap structure 
obtained at T = T c is fairly robust as a function of T 
and, in particular, if nodes exist near T c they can sur- 
vive as T —> 0; ii) at the zero-points A' = Ao, there are 
gapless quasiparticle excitations with unusual anisotropic 
dispersion. This leads to a distinctive power-law temper- 
ature behavior of many thermodynamic quantities; iii) At 
T = these zero-points form a surface of quantum phase 
transitions with novel and peculiar anisotropic scaling. 
Consequently, the familiar Simon-Lee scaling 2 ^ of Dirac- 
like nodes is superseded by a different form of scaling - 
this provides us a with a new diagnostic tool to apply to 
the phenomenology of pnictides. 



A e at (j>f,0), (0,p F ), (-p F ,0) and (0,-p F ) (using coor- 
dinate system, rotated by ir/4 with respect to BZ axes, 
and where pp is the Fermi momentum for the electron 
band). The same structure is reproduced at other corners 
of the BZ. It is important to stress here that, while ([1]) 
is a simplified phenomenological model, the physics and 
the results that are the focus of this paper will remain 
unchanged when a more realistic description is employed. 




II. THE STRUCTURE OF THE GAP 
FUNCTION 

The two-band model reflects the basic features of iron- 
pnictides: a hole-like and an electron-like band at the 
center and at the corner of a square Brillouin zone (BZ), 
respectively. These bands are coupled via interband pair- 
ing interaction N(0)V(k, k') = A +A„ cos 48, where N(0) 
is the density of states (DOS) at the FS, Ao and A ra 
are coupling constants and 9 is an angle on the elec- 
tron FS (using coordinate system centered at the elec- 
tron pocket). Defining gap functions Ah (6) = Ah and 
A e (0) = A + A' cos 49 and introducing intraband repul- 
sion fj, results in the gap equations: 
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where A is a high-energy cut-off and E e = 
yj^ 2 + A 2 ; Eh = \/£, 2 + A\. Since all interactions are 
repulsive (A e ) FS'Ah < is a necessary condition for non- 
trivial solutions of Eqs. ([T]) to exist ((•••) fs is the an- 
gular average over the electron FS). Such sign-switching 
state is a likely prospect for pnictides 2 . For T — > T c , these 
equations can be linearized with the following result (see 
also Ref. 1231 ): when A' = Aq, there are four zeroes in 



FIG. 1. A h (T)/T c as obtained from solving eqs. (TjJ. The 
solid red line is for parameter values \i — 0.05; Ao = 0.2 and 
A n = 0, and corresponds to the BCS result for a single band. 
Increasing \ n = from zero leads to decrease in Ah/T c - 
shown for A n = 0.1, A n = 0.15, A n = 0.2 and A„ = 0.3 
(the blue component increases with increasing An). Note the 
negative sign of - consequence of choosing Ao positive. 




FIG. 2. Ao(T) obtained from solving eqs. |T}. The parame- 
ters, color coding and normalization are the same as in Fig. 
□ 

The first such result follows from solving eqs. ([T]) at 
arbitrary T, sl complicated task which we accomplish nu- 
merically. Representative results for Ah, Ao and A' and 
are displayed in Figs. Q] [2] and [3l For A„ = the model 
reduces to the case of two uniform gaps which follow 
(for equivalent bands) the BCS temperature dependence 
(a known result for multiband superconductors 2 ^^). A 
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FIG. 3. A' as function of temperature. The parameters, color 
coding and normalization are the same as in Fig. [1] 



nonzero value of A n generates a finite A', and pushes 
A /T c and A'/T c below their BCS values (of course, 
in absolute units T c ,Aq and A^ all go up). Evidently, 
the temperature dependence of different gaps is similar. 
Most importantly, this implies that A'/A has a rather 
weak T dependence and a nodal state dXT — T c will thus 
likely remain nodal as T goes to (shown on Fig. 0}. 
Remarkably, A'(T)/Aq(T) can both increase or decrease 
or can even exhibit a non-monotonic behavior, testifying 
to the strong non-linearity of the gap equation. Never- 
theless, the overall changes are generically quite small. 
Thus, the T = T C phase diagram survives at T — with- 
out major changes (Fig. [5]shows our phase diagram, with 
A fixed). We also point out that A (T -> 0)/A h (T -> 0) 
is not equal to A (T T c )/A h (T -> T c ) as expected for 
isotropic two-band superconductor—. 

The above results are surprising since the nodes are 
only accidental and therefore sensitive to disruption by 
the strongly non-linear nature of the gap equations ([1} . 
Their robustness can be viewed as a posteriori justifi- 
cation for the use of nodal or near nodal states in low 
T calculations, since previously these states were self- 
consistently obtained only at T c . Note also that the three 
coupling constants of our model lead to a certain am- 
biguity: a nodal state can be created by large A n or 
by strong fjb, since intraband repulsion suppresses the 
uniform gap stronger than the oscillating one. These 
two distinct mechanisms of generating nodes - a very 
anisotropic interband interaction and a strong intraband 
repulsion - are both present in pnictides, and reflect gen- 
uine physics of these multiband superconductors. Fur- 
thermore, because the number of gaps in real materials 
is (at least) four, there is additional frustration associated 
with interband pair scattering between the gaps with the 
same sign. This frustration effect also favors the nodal 
components 22 , and should be included in a more realistic 
four or five-band calculations. 

The second key result follows from the first: consider 
the surface A' = A in the (A„,/z, T) space (with Aq 
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FIG. 4. Ratio A'(T)/A'(T), arising from solution of Eqs. Q 
for n = 0.05; A = 0.2 and A n = 0.1 , 0.15, 0.2, 0.3, respectively. 
This ratio does not have a simple temperature behavior, but 
the changes are generically quite small. This means that the 
nodes, present at T — T c , likely survive as T — > 0. 



fixed), and ignore A^ since its contribution is exponen- 
tially suppressed at low T . On this surface the gap in the 
electron band develops four zeroes, which, for A' > Ao, 
split into eight nodes. At T = 0, this surface defines 
a line of quantum phase transitions (see Fig. [5]), along 
which the low energy quasiparticle spectrum suffers a 
non-analytic transformation. For A' < Aq, the quasi- 
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FIG. 5. The A'/Ao isolines for various {i and A„, with Ao = 
0.2. All of them converge at = Ao point at the j/-axis. 
The color axis represents A'/Ao. For numerical convenience 
T is kept finite, but extremely small (T = I0 -5 T c ), and is 
practically zero. Inset: color map of T c (in arbitrary units) 
for the same values of /i and A n . 



particle excitations are gapped whereas, for A' > Aq, 
one enters the nodal region, with gapless Bogoliubov- 
deGennes (BdG) quasiparticles described by a Dirac-like 
Hamiltonian 27 . At the critical line A' = A separat- 
ing the two regimes - where the pairs of nodes merge 
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into zeroes - the quasiparticle dispersion has a peculiar 
character: linear along the Fermi velocity and quadratic 
perpendicular to it (see Fig. [6|). The resulting spectrum 
gives rise to an unusual critical scaling of low T thermo- 
dynamics, which dominates the phase diagram surround- 
ing the critical line. 



III. ZEROES IN THE GAP FUNCTION AND 
THEIR SIGNATURES IN QUASIPARTICLE 
PHENOMENOLOGY 

To explore this peculiar criticality at the A' = Ao 
line, we expand the gap function in the BdG Hamiltonian 
around one of the zeroes, say (j?f,0): 
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Eq. 



^ results in an anisotropic dispersion E 
for 



\J (vfPx) 2 + (8A Py/p|0 2 f° r quasiparticle energies. 
This leads to anisotropic scaling and different low-energy 
(E <C A Q ) scaling lengths along x and y directions: vpjE 
and y/ vpa/E, respectively, where a = (8Aq/vfPf)- In 
the Lagrangian nomenclature, the scaling dimensions at 
T = are [y] = —1, [x] = —z x , [r] = — z T , where 

~ ~ o31 



We have calculated the total DOS - N tot (E,T) = 
N e (E, T) + N h (E, T) - by using the gap functions shown 
in Figs. [[J M and EI The results (at T = 0) are dis- 
played in Fig. [7J Clearly, the DOS exhibits a complex 
behavior with several prominent features. First, the two 
coinciding (at A n = 0) BCS singularities in N e and Nh 
at E = A/j = A e , split for a nonzero A„. In addition, 
a new peak appears in N e , due to the introduction of a 
new energy scale A'. The positions of the three peaks 
evolve with increasing A„. For A' > A , the low-energy 
excitations exist all the way down to E = 0, as expected 
for nodal order parameter. 
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FIG. 6. The quasiparticle dispersion near accidental zeroes 
and nodes. For strongly anisotropic gap, |Ao — A'| <C Ao, the 
"zero-point" critical scaling dominates the nodal contribution. 

We can now evaluate the DOS of BdG quasiparticles: 

E 



N e (E) = Re 



VE (3) 



FS 



which is different from the familiar nodal DOS Na(E) ~ 
E. This increase of DOS relative to the nodal case is 
caused by the parabolic part in E. An exact expression, 
in terms of complete elliptic integrals^, is: 



N ^ = ^\li K ^ where W 1 + ^o ;(4) 

where N z is the number of gap zeroes, and as expected, 
for E -> 0, N e (E) -> \/E. 



FIG. 7. The total DOS, calculated using the gap functions for 
X„ — 0.1,0.15,0.2,0.3, respectively. In each panel the BCS 
DOS (An = 0) is included for comparison. 

Having determined DOS, we are now in position to 
extract the low T dependence of various thermodynamic 
quantities. The specific heat can be written as: 



C(T) 



El 



dE- 



E 2 N(E) 
cosh 2 (E/2T) 



N Z T 3/2 . (5) 



Similarly, we can obtain the temperature dependence of 
the spin susceptibility and the superfluid density: 



Xs(T)/Xn = 1- Ps{T)/p 



dE- 
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from which we can obtain the leading temperature- 
dependent term of the penetration depth: 



A(T) 



^ 47re 2 p s (T) 
The nuclear spin relaxation rate is: 

(T\)„ ^N, [°° dE _ N{Ef 



A(0) + AiT 1 / 2 
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For the thermal conductivity we need the quasiparticle 
scattering rate, which can be written as: 



1 



AqAi, 



where t~ an d iV n are the normal state scattering rate 
and DOS respectively. The thermal conductivity in the 
superconducting state is: 
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where ki is the ith component of the quasiparticle mo- 
mentum. Integrating around the zeroes and keeping only 
the term with the lowest power of the temperature, we 
get for the thermal conductivity, normalized by the nor- 
mal state value: 



(«n/T) 



J^2 



dEE 5 ' 2 
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N Z T. (9) 



All of the derived scaling forms strictly apply only at 
the critical line. Of course, tuning precisely to a point 
on this line is not an easy task experimentally, and this 
means that the asymptotic low temperature behavior of 
the thermodynamic quantities will ultimately be gov- 
erned either by the nodal Dirac or by thermally excited 
gapped BdG quasiparticles, depending on which side of 
the line the system finds itself. Importantly, however, as 
long as one is in the vicinity of the quantum critical line, 
this asymptotic behavior will be restricted to a very nar- 
row low temperature range, T <C | Ao — A' | , and there will 
be a large wedge-shaped crossover regime, dominated by 
the quantum phase transition, for which the temperature 
scalings are governed by Eqs. ([S])-®. 

We supply an illustration of the above ideas by calcu- 
lating the spin susceptibility as A' — > Ao, and the result 
is shown in Fig. [3] The low-temperature part of Xs(T) is 
clearly enhanced, and close to A' = A (for A„ = 0.15) 
the T 1 / 2 crossover behavior is expected. 



IV. GAP ZEROES IN FINITE MAGNETIC 
FIELD 

Next, we turn to the problem of a superconductor with 
gap zeroes in an external magnetic field H. Finite H 
modifies the critical line in the (A n ,/x) plane, but, for 
H far below the upper critical field H C 2(0), such modifi- 
cations are negligible. The presence of gapless fermions 
makes the finite H problem rather nontrivial, due to the 
singular scattering of fermions from vortices 3 - 3 -. Useful 
results, however, can be extracted from general scaling 
arguments; in nodal d-wave superconductors such argu- 
ments are behind the well-known Simon-Lee scaling 2 ^. 

We concentrate on intermediate fields (H c i <C H <C 
Has), where the vortex spacing is large and we can 




FIG. 8. The spin susceptibility, calculated using the temper- 
ature dependence of the gap functions for A n = 0,0.1, and 
0.15. The enhancement of the low-temperature part is clearly 
visible. The case A n = 0.15 lies almost on the critical line and 
thus Xs(T) ~ T 1 / 2 is expected as a crossover behavior at low 
temperatures. 



assume a uniform field and a constant gap ampli- 
tude. Near the (p F ,0) zero, the continuum limit of 
the gauge- invariant gap function takes the form A(r) = 
(8/p£O{0„, {d v , A(r)}} + (2i/p%)A(r)d$<j>, where A(r) = 
|A|e^ < - r ' ) and {,} signifies full symmetrization 3 - 3 -. This 
allows us to write the low-energy BdG Hamiltonian as: 
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(10) 



We now choose Landau gauge A x — Hy and rewrite 
vf(Px + ifi/c)Hy) in dimensionless form by introduc- 
ing a length d and defining variables x' = x/d and 
y 1 = yd/l 2 , where I = ^/c/eH is the magnetic 
length. The off-diagonal terms in "HedG a- re similarly 
rewritten as (8A /p 2 F )[{d y , {d y , e*+}} + (i/4)e^(d 2 ^)} = 
(w F ad 2 /Z 4 )[{5y,{9^,e^}} + (i/4)e^(5 2 ,0)] and its com- 
plex conjugate. The choice d = (^/a) 1 / 3 allows us to 
reexpress the quasiparticle energy spectrum as {E n } — > 
(v F /d){e n } = v F (a/l^/*{E n } = v F (ae 2 H 2 /c 2 ) 1/3 {e n }, 
where {s n } are numbers which in principle can depend 
only on H/(c/ea 2 ) = H/H a <C 1. If they do, such 
"anomalous" scaling would not be surprising in our case 
of singular anisotropy. Adding the contributions from all 
zeroes restores the x O y symmetry. 

Using the above rescaling, we write DOS N(E, H) as: 



U^(ae 2 H 2 /c 2 ) 1 ^( 



E _ H 

v F {ae 2 H 2 /c 2 ) 1 /^ ~H~J ^ ' 

where Af(u; w) is the scaling function. If H induces fi- 
nite DOS at low E, we can follow earlier discussion to 
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obtain N(E -+ 0,H) ~ H 1 / 3 ^', where T) > -1/3 is 
the "anomalous" dimension, defined by Af(u = 0;w — > 
0) = const, x w n . r) is clearly a possibility given 
the fact that the vortex lattice structure in exp(i</>(r)) is 
determined by all electrons and does not conform to the 
anomalous scaling of individual zero-points. The precise 
value of r\ follows from direct computation^. 

The critical "zero-point" scaling also works when both 
H and T are finite. The internal energy density can be 
written as 

U(H,T) = ^£„(fl)/(£„(ff)/T) ~ (12) 

n 

~ N z H 5 / 3 T u (T/(v 3 F ae 2 c- 2 Hy/ 3 ;H/H a ) , 

where / is the Fermi function and J- u is a scaling function. 
Here we converted ^ n — > (V/Att 2 ) J dp x dp y , dp x dp y ~ 
(\/l 2 )dp x idp y £^. From (|13p. the specific heat per unit 
volume is: 

C(H,T) ~ N z H 5 / 3 drMT/(v 3 F ae 2 c- 2 H 2 y/ 3 ;H/H a ) ~ 
- N.HF.iT/ivla^c^H^-H/Ha) , (13) 

where J" C (T '/ {v F ae 2 c~ 2 H 2 ) 1 / 3 ] H/H a ) is the correspond- 
ing scaling function. 

In general, evaluating ,F„ or J-" c is a difficult task. How- 
ever, the needed limits of these functions are readily de- 
duced: Eq. (|SJ| mandates .F c (a — > oo; tu — s- 0) — > a 3 / 2 , 
where a = T /{v F ae 2 c~ 2 H 2 ) 1 / 3 ' . In the opposite case, 
a — > 0, we can use the constant low E DOS, induced by 
H, and observe that © gives C(H,T) ~ TH 1 ^, and 
thus J-" c (a — ► 0;w <C 1) — > am 11 . Note that these scal- 
ing arguments apply only to the critical BdG fermions 
inhabiting the zeroes. We are not including the con- 
tributions from other parts of the system, like possible 
localized states in the vortex cores - this contribution is 
subleading at low E, T and H. 

Again, we want to emphasize that although precise 
tuning to a state with gap zeroes seems very unlikely, 



close to such state and at not too small T and H the 
scaling associated with gap zeroes dominates. Thus, in 
a state with accidental nodes the Simon-Lee scaling can 
be completely unobservable. 



V. CONCLUSIONS 

In summary, we have studied an anisotropic two-band 
model of iron-based superconductors. By explicit solu- 
tion of the gap equations, we uncovered the robust na- 
ture of the accidental nodes and zeroes in the gap func- 
tion. We have discussed in some detail the quantum crit- 
ical line where the gap zeroes first appear. It projects 
a considerable influence over the phase diagram and is 
characterized by a peculiar form of anisotropic critical 
scaling, qualitatively distinct from the familiar Dirac- 
Simon-Lee scaling. Irrespective of whether one can tune 
in precisely to this quantum critical line in a particular 
experiment, as long as the gap is strongly anisotropic 
(|Aq — A' | <C A ), the "zero-point" critical scaling will 
dominate the quasiparticle thermodynamics and trans- 
port over a wide region in a phase diagram, overwhelming 
the contribution from accidental nodes. 
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